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Abstract 

In the semi-classical approach to the Skyrme model, nuclei are approximated by 
quantum mechanical states on a finite-dimensional space of field configurations; in 
zero-mode quantization this space is generated by rotations and isorotations. Here, 
simulated annealing is used to find the axially symmetric Skyrme configuration which 
extremizes the zero-mode quantized energy for the nucleon. 

1 Introduction 

The Skyrme model is an effective theory of pions and nucleons. It is a non-linear field theory 
in which nuclei correspond, classically, to topological soliton solutions, called Skyrmions. 
The model is non-renormalizable and so the approach usually taken is to reduce to a finite- 
dimensional space of Skyrme configurations before quantizing. In this approach, quantum- 
mechanical states on a space of topological charge B Skyrme configurations model baryon 
number B nuclei at low energies. 

The Skyrme Lagrangian is acted on invariantly by rotations and isorotations; this 
generates from a single Skyrme configuration a space of energy-degenerate configurations. 
This space is known as the space of zero modes and quantization on this space is called 
zero-mode quantization. This approach began with the seminal paper £Q where the nucleon 
and delta are constructed as quantum mechanical states on the zero-mode space generated 
from the classical minimum energy configuration of unit baryon number. 

In this paper, the spin-half, isospin-half quantum Hamiltonian is calculated on the 
zero-mode space generated from a general axially symmetric configuration and simulated 
annealing is used to find the Skyrme configuration that minimizes the energy of the lowest 
state. 

This approach has previously been considered for the zero-mode space of a general 
spherical symmetric configuration Ej • In this space the quantum Hamiltonian reduces 
to a scalar, here it is a matrix. Aspects of our approach are also shared with a very recent 
paper JI| which appeared when our paper was in preparation. The quantum Hamiltonian 
used in that paper is a scalar ansatz motivated by the spherical Hamiltonian used in pQ 
and differs from the Hamiltonian we derive here. In fact, we will see that this will not 
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make a significant difference, the numerical results obtained for the nucleon in both papers 
do not differ much from each other or from what would be calculated using the classical 
minimum. However, the approach here is more direct than the approach described in 
and can be generalized to higher charge nuclei. 



2 Quantization procedure 



Written in terms of the vector currents = d^UW of an SU(2) field U (x), the Skyrme 
model has the Lagrangian 
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-^Tr (R^) + ^Tr ([R„ R u ] [R», it"]) + Tr (U - 1) 
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where m n , F n and e are parameters that are adjusted to fit experimental data. Using F^/Ae 
as our unit of energy and 2/eF T as our unit of length we obtain 
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The Skyrmion mass for a static field U s (x) can be derived from this Lagrangian and is 
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We wish to quantize the rotational and isorotational degrees of freedom. Rather than 
acting on a specific Skyrme configuration, we want to consider the zero-mode space of fields 
generated from a general static configuration U s (x) by isorotation C and rotation D: 



U(x) = CU s (x D ) c\ 
where C is in the 2x2 representation, 
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and Dij is a three-dimensional matrix representation of D. Since rotation and isorotation 
are symmetries of the original Lagrangian, these configurations are all energy-degenerate. 
The effective Lagrangian on this restricted space of configurations can be calculated by 
allowing C and D to depend on time, giving L = —M + L rot , where L rot is the kinetic 
Lagrangian 

l l 
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L ro t — —QiUijQj + -^iVijUj — UiWijQj, 

with the rotational and isorotational angular velocities uj and f2 given by 



iTr ( (TiD^D 
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and the moment of inertia tensors and H^,- given by the following integrals 
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Tr(T t T J ) + -Tr([R k ,T l }[R k ,T J }) 
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Tr (T,i? m ) + - Tr ([R k , T t ] [R k , R m ]) 



and 
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where the cr; are the usual Pauli matrices. 

In this paper, we will restrict our discussion to axially symmetric Skyrmion solutions. 
Numerical simulations indicate that there is also a reflection symmetry in the xy-plane 
and so the principal axes of inertia can be taken as the standard orthogonal axes, with 



Uij — Vij 



W, 



where i ^ j, and we can set Uu = Ui and so forth. By inspection, 



cylindrical symmetry in the xj/-plane will also mean 

U x = U 2 , Vt = V 2 , Wt = w 2 . 



(10) 



We can use axial symmetry to establish an additional identification between the normal 
moments of inertia (see Appendix): 



u 3 = v 3 = w 3 

Applying these restrictions to L rot © we get 
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Lrot = ~ (fi? + nl) U 2 + i (uj + u\) V 2 + ^ (Q 3 - oo 3 f U 3 - (Qtoj, + n 2 oo 2 ) W 2 (12) 

or, written as a sum of complete squares, 
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The rotation and isorotation angular momentum vectors L and K canonically conjugate 
to to and f2 are 
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K 



dL 



rot 



dL rot 

on 



(v 2 u x - w 2 n l7 v 2 u 2 - w 2 n 2 , u 3 (u 3 - n 3 )) , 
(u 2 Qi - w 2 ui, u 2 n 2 - w 2 u 2 , -u 3 (u 3 - n 3 )) . 
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Note that L 3 =—K 3 : upon quantization, this is the axially symmetric condition expressed 
in operator form. Our approach will be to find the minimum energy Skyrmion such that its 
energy eigenstate is also a null eigenstate of L 3 + K 3 . From (|T3|) and (|14)1. the Hamiltonian 
for the rotational and isorotational degrees of freedom can now be calculated: 
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rot 
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For a spin-n,isospin-n particle, the angular momentum operators L can be written as 
(2n + l) x (2n+l) dimensional matrix representation E^, E^ , E3 of SU(2), as can the 
isospin operators Ef , Ef , E3 of K for an isospin-n particle. If we define our quantum 
state using the \l, I3) <S> \k, k 3 ) basis, I, I3 and k, k 3 being the quantum numbers for L and 
K respectively, we can embed L and K into the resulting SO(3) L x SO(3)^ direct product 
space: 



L \-> KZ L ®I 2n+1 



(16) 



where I 2n +i is the (2n + l)-dimensional identity matrix. 

We can now find the lowest energy nucleon state; first we insert the spin-half, isospin- 
half matrix representation of L and K into the Hamiltonian Ijl5j) to get: 
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where 



and 
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There are two eigenvectors of H which are also eigenvectors of L 3 + K 3 with eigenvalue 
zero: 



[0, 0> 
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1. ,11 ,11 ,1 l.\ K 1 , 

1,0) = -= I-,--}® |-,-) + |-, -}® |-,--) , Si,o = t («i + «2) (21) 



and 

v^V' 2 ' 2' ~ '2' 2' ' '2' 2' ^ '2' 2^ ' ~ X,U ~ 4 

The first eigenvector is a spherically symmetric state and also has the lower energy since 
K\ and k 2 are always positive; hence 
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and is therefore the energy of an axially symmetric nucleon. We see that in the spherically 
symmetric case U 2 = U 3 = W 2 = V 2 , and E N reduces to the rotational energy formula 
obtained in £Q: 

Br .j| ((i+1H y- (23 ) 

where A = | {U\ + U 2 + U3) = U3. The energies of all other eigenstates go to infinity in 
the spherically symmetric limit. 

The energy, M+E N , of the quantum state would be difficult to extremize using gradient- 
based methods; instead, simulated annealing [S] is used to find the Skyrmion configuration 
that minimizes this energy. Since the configuration is assumed to be axially symmetric, 
the cartoon method t6. is used. The configuration is annealed on a quarter-disk two- 
dimensional lattice with a radius of 250 lattice points and a lattice spacing of 0.06. A 
variant of the Adaptive Simulated Annealing probability distribution |7j is used for the field 
perturbations; this seems to improve the speed of convergence and allows an exponential 
cooling schedule. The algorithm needs an initial configuration to perturb; any configuration 
of unit baryon number will suffice, and the ansatz given in [B] is probably the easiest to 
implement. 

In Fig.l and Fig. 2, the results of our simulations are shown in comparison with results 
obtained using the rigid body approximation used in 0. The pion mass was set to its 
experimental value of 138 MeV in Fig.l, whereas in Fig. 2 it was set to the larger value of 
345MeV suggested in j^; we see that the nucleon deformation only becomes noticeable at 
high values of the pion mass. 



3 Discussion 

Although it has appealing mathematical and physical properties, the Skyrme model has 
only been modestly successful in modeling nuclear dynamics. Perhaps one of the biggest 
challenges is that the zero-mode quantization of the classical minimum fails to give even the 
correct lowest energy state for many values of the baryon number [TTH ITT] . It is possible 
that this is because inappropriate parameter values are used to calculate the classical 
minima: for example, in it is suggested that using a non-zero, or even an unphysically 
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Figure 1: Plot of F^, as a function of e, for which the energy M + En is equal to the 
nucleon mass 939 MeV; the pion mass parameter m n set to its experimental value of 138 
MeV. Our results (bold circles) are compared with those obtained using the rigid body 
approach taken in (solid line). 




Figure 2: As Fig. 1, but with the pion mass parameter set to the value of 345 MeV 
suggested in [3]. Our results (bold circles) are compared with those obtained using the 
rigid body approach taken in [0] (solid line). 



6 



large, pion mass will significantly affect the structure of higher charge Skyrmions and alter 
the symmetries on which the quantization is based [T3]. Another possibility is that the 
classical minimum is inappropriate and that the energy minimum of the effective quantum 
Hamiltonian should be used for zero-mode quantization. In this paper this has been done 
in the simple case of a axial symmetric nucleon; the higher charge cases will be more 
challenging computationally but should follow in a similar way. 
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Appendix: Identification of normal moments of inertia 



Following [Mj, we establish a relation between U 3 , V 3 and W3 for axially symmetric fields. 
We first express the quantity e 3 jkXjRk in polar coordinates: 

e 3jk XjR k = (x x V) 3 UU ] 

- %«■ <*> 

We can then identify e^x^Ri and — | [03, U] by looking at the general form of an axially 
symmetric field with unit baryon number: 

— i<7rtd> . r / \ i<j*)d) 

U = e ^% 2/(r '* )niCTl e^ (25) 
and taking its derivative with respect to </>: 

dU_ = _w 3 ^_^ ei/(riZ)e ^j + ^z^± elf{r:Z)e ^ ^3 

= ( 26 ) 

We see the expressions for the inertias C/3, V3 and W3 in the moment of inertia integrals 
(JHJ) differ only in the terms e^kXjRk and — | [a^, U] W, and so 

U 3 = V 3 = W 3 . (27) 
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